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Abstract. We discuss trivial deformations of the canonical Poisson brackets associated 
with the Toda lattices, relativistic Toda lattices, Henon-Heiles, rational Calogero-Moser 
and Ruijsenaars-Schneider systems and apply one of these deformations to construct a new 
trivial family of noncommutative integrable systems. 
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1 Introduction 

Let us consider some smooth manifold M with coordinates x%, . . . , x m and a dynamical system 
defined by the following equations of motion 

Xi = Xi, i = l,...,m. 

We can identify this system of ODE's with the vector field 



which is a linear operator on a space of the smooth functions on M that encodes the infinitesimal 
evolution of any quantity 



In Hamiltonian mechanics one of the fundamental axiom is what we can call the energy paradigm 
that can be stated as follows: "For every mechanical system there is a function defined on its 
space of states, called mechanical energy or Hamiltonian H of the system, containing all its 
dynamical information" . 

According to this paradigm any function H on M generates vector field X describing a dy- 
namical system 



Here dH is a differential of H, and P is a bivector on the phase space M. By adding some 
other assumptions we can prove that P is a Poisson bivector. In fact, it is enough to add energy 
conservation 





X = X H = PdH. 



H = X H {H) = (PdH,dH) = 
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and compatibility of dynamical evolutions associated with two functions H\ % 

X Hl (X H2 (F)) = X H2 (X Hl (F)) + X XHi(H2) (F), 

see [17, 20] and references therein. 

In bi-Hamiltonian mechanics [25] we are looking for another decomposition of the given vector 
field X, 

X = PdH = ftP'dHx + ■■■ + f m P'dH m , 
by commuting Hamiltonian vector fields 
Xk = P'dHk 

generated by integrals of motion Hi, . . . , H m and some Poisson bivector P' compatible with P. 
These Poisson bivectors P' can be divided into two groups of trivial and nontrivial deformations 
of canonical Poisson bivector P, see [3, 9, 23, 25, 33] and references therein. 

Supposing that P' is a trivial deformation of canonical Poisson bivector P completely defined 
by Hamiltonian H we can join the geometry of the phase space M and the energy paradigm. In 
this case, second Poisson bivector 

P' = C Y P (1.1) 

is a Lie derivative of the canonical bivector P along the Liouville vector field 

Y = AdH, (1.2) 

where A is a 2-tensor field acting on the differential of the Hamiltonian. Remind, that the 
Lie derivative P' (1.1) is a trivial deformation because it is 2-coboundary and simultaneously 2- 
cocycle in the Poisson-Lichnerowicz cohomology defined by canonical Poisson bivector P [9, 23]. 

The main aim of this note is to show 2-tensor fields A associated with some well-known 
integrable systems and to prove that these tensor fields may be useful to construct new integrable 
systems. For example, we discuss a new trivial family of integrable noncommutative three 
dimensional systems, which includes deformations of the rational Calogero-Moser system with 
three particle interaction. 



2 Integrable systems on cotangent bundles 

Let us consider canonical Poisson bivector on the symplectic manifold M = T*Q 

(2 - 1) 

which is the one mostly used in Hamiltonian mechanics [1, 2]. Here q% are local coordinates 
describing a point q on a smooth manifold Q, and pi, the canonical conjugate momenta, are 
local coordinates describing covectors on such manifold, i.e. points p in the cotangent bundle 
T*Q of Q. The corresponding Poisson bracket looks like 

{Pi,<lj} = 8ij, {q%,Qj} = {PhPj} = 0- ( 2 - 2 ) 
In local coordinates x on M the Lie derivative of a bivector P along a vector field Y reads as 

3 ^ V dx k dx k oxkj 
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and the Schouten bracket [A, B] of two bivectors A and B is a trivector with the following entries 

dim M , qj^ x 

[A,B] ijk = - ^ (B mk ^+A mk ^ + cjcle(i,j,k)j. (2.3) 



In our case M = T*Q, dimM = 2n and x = (q,p). 

Let us consider natural Hamilton functions on M = T*Q 

H = T(q,p) + V(q), (2.4) 

which are the sum of the geodesic Hamiltonian T and potential energy V(q). According to [15, 
24, 31, 32], for natural Hamiltonians there is other representation for the vector field Y (1.2) 



Y 



A 0\ fdq 
Uj \dp 



Here we use matrix notation of tensor objects in which, for instance, canonical Poisson bivec- 
tor P (2.1) looks like 

*=(-°* »)• M 

Here E is a unit matrix. 

If n = and A is the conformal Killing tensor of gradient type or Yano-Killing tensor on Q, 
one gets second Poisson bivector P' = LyP associated with Hamilton functions separable in 
orthogonal coordinate systems on Q. In this case eigenvalues of the the Nijenhuis operator 

N = P'P- 1 , (2.6) 

which is also called the hereditary or recursion operator, are variables of separation. In order 
to get integrals of motion H k from N we have to extend the initial phase space [18] or to fix 
separated relations [31]. 

The new idea is that we can substitute the arbitrary Hamilton function H (2.4) and the 
2-tensor field A into the definition (1.1) and try to find the Poisson bivectors P' solving the 
equation 

[P',P'] = [C AdH P,C AdH P] = 0, (2.7) 

where [•, •] is the Schouten bracket defined by (2.3). In this case P' will be the Poisson bivector 
compatible with P and we will say that M = T*Q is the bi-Hamiltonian manifold [25]. The 
next step is a search of integrals of motion for this Hamilton function H. 

If the recursion operator iV at every point has n distinct functionally independent eigenvalues, 
we can say that M is a regular bi-Hamiltonian manifold. If the recursion operator N does not 
have this property then we can say that bi-Hamiltonian manifold M is irregular [24, 31]. So, 
there are three different cases: 

1) recursion operator produces the necessary number of integrals of motion; 

2) recursion operator generates variables of separation instead of integrals of motion; 

3) recursion operator produces only part of the integrals of motion or variables of separation. 

In the third case we have to complement the recursion operator with some additional information 
in order to get integrals of motion. Namely this property allows us to get noncommutative 
integrable systems, which will be considered in Section 3. 

Now let us show a collection of tensor fields A associated with some well-known integrable 
systems. 
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2.1 Toda lattice 

Let us consider the following tensor field A depending only on q variables 



A 



B + 2D- B 







P, 



B-B 1 
where B is a strictly upper diagonal matrix 



B 



/0 1 1 
1 



... ix 
... i 

1 
0/ 



e ^ ' 



(2. 



and Z) is a diagonal matrix 

n 

D = diag(gi,g 2 , .■■,q n ) = ^%e. 



i=l 



Matrices eij are n x n with only one non zero (ij) entry, which equals to unit. 
Substituting this tensor field A and a natural Hamilton function 



i=l 



into the definition of P' one gets a system of equations (2.7) on V(q). One of the partial solutions 
of this system is the Hamilton function for the open Toda lattice associated with A n root system 



n-l 



a G 



i=l i=l 

Traces of powers of the corresponding recursion operator N (2.6) 
H k = trN k , k = l,...,n, 



(2.9) 



are functionally independent constants of motion in bi-involution with respect to both Poisson 
brackets 

{H i ,H j } = {H i ,H j }' = 0. 

This Poisson bivector P' was found by Das, Okubo and Fernandes [8, 10]. 
In generic case we can use a more complicated tensor field 



A 



B + D 







C, 



P, 



where entries of D are linear on qi and B and C are numerical matrices. Here B is an arbitrary 
matrix, whereas A and B satisfy to algebraic equations which may be obtained from (2.7) at 
V(q)=0. 

Using this tensor field A we can get the recursion operators which produce either integrals of 
motion for the periodic Toda lattice [14] or variables of separation for the Toda lattice [34]. In 
similar manner we can consider the Toda lattices associated with other classical root systems [31]. 



On a Trivial Family of Noncommutative Integrable Systems 



5 



2.2 Relativistic Toda lattice 

If we substitute the Hamilton function H(q,p) and the following tensor field A 

{-E Oj \0 -E 

where E is a unit matrix and B is given by (2.8), into the definition of P' (1.1) we will obtain 
a system of equations on H. One of the solutions is the Hamiltonian of the open discrete Toda 
lattice associated with A n root system 

n 

H = J2i c i + d i ), (2.10) 

where and d{ are the so-called Suris variables 

Ci = exp(pi - qi + q i+ i), d» = exp(pi), q = -oo, q n+1 = +oo. 

Traces of powers of the corresponding recursion operator N (2.9) are integrals of motion in 
bi-involution with respect to both Poisson brackets. Namely this Poisson bivector P' (1.1) is 
discussed in [26, 29]. 

Remind, that according to [29] there is an equivalence between the relativistic Toda lattice 
and the discrete time Toda lattice. Namely, substituting 



v =0 I 1 ln / < 1 + ex P(gj~gj-i) > \ 
3 3 2 \1 + ex V (q j+1 - gj ) J 



in (2.10) one gets standard Hamiltonian for the relativistic Toda lattice 

n ~ l r 1 V 2 

H = J2 ex P(0j) [l + exp(gj - + exp(g i+ i - q^)] 

Transformation [0j ■ qj) —> (pj ■ qj) is a canonical transformation. 
As above, two numerical matrices B and C in the tensor field 



A 



B 
C 



allow us to get recursion operators N = P'P 1 which generate either integrals of motion for the 
periodic relativistic Toda lattice or variables of separation [22]. 

2.3 Henon— Heiles system 

At n = 2 we can introduce the following linear in momenta tensor field A 



where 



B = ( 2 Qm QiP2\ c _ ( fi(q)pi + f2{q)P2 



\qiP2 Q2P2J ' V h(q)Pl + fa(q)P2 
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Substituting this tensor field A and a natural Hamilton function 
H^pl+pl + Viq) 

into the definition of P' one gets a system of equations (2.7) on V(q) and functions f k (q)- The 
resulting system of PDE's has two partial polynomial solutions 



V(q) = c iq2 {Sqf + IQql) + c 2 (2ql + + C3</2, c k £ M, 



and 



V(q) = cx {qf + 6qjq 2 + 8q 4 2 ) + c 2 {<£ + Aq 2 ) + 4. 

Second integrals of motion H 2 = tr N 2 are fourth order polynomials in momenta. 

So, one gets the Henon-Heiles potential and the fourth order potential [12] as particular 
polynomial solutions of the equations (2.7) associated with tensor field (2.11). 

Using slightly deformed tensor field A we can get the same systems with singular terms [15] 
and their three-dimensional counterparts [31]. 

2.4 Rational Calogero— Moser model 

Following [24] let us consider tensor field A, which is proportional to P 

A = p(q,p)P, 
where p(q,p) is a function on M. If 

A = (piqi H h p n q n )P, P = PiQi H h PnQn, (2.12) 

then equations (2.7) have the following partial solution 
1 n 2 n 1 

where g is a coupling constant. It is the Hamilton function of the n-particle rational Calogero- 
Moser model associated with the root system A n . 

The corresponding recursion operator ./V (2.6) generates only a Hamilton function 

trN k = 2H k , k = l,...,n, 

that allows us to identify our phase space M = R 2n with the irregular bi-Hamiltonian manifold. 
In this case [24, 31] integrals of motion are polynomial solutions of the equations 

PdH = -\p'd In H k , k = l,...,n, (2.14) 
k 

which have two functionally independent solutions for any k > 2. It is easy to see that the 
functions 

TT-l/m 

C km = (2.15) 
are Casimir functions of P', i.e. P 'dC km = 0. 
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Some solutions of equations (2.14) coincide with the well-known integrals of motion 

m = 1, . . . , n — 1, 



~~ m) 1 ^ " ' 1 ^ I I ' 

v v ' 

m times 

obtained from the conserved quantity 

( g 2 1 <9 2 \ - 

Jn s exp l -t g te-^) 2 ^- ) n 



by taking its successive Poisson brackets with ^ [11]. These n solutions, including J2 = H, 



i=l 



are in involution with respect to the Poisson brackets (2.2). 
Other n — 1 functionally independent solutions of (2.14), 



K m = mgiJm - g m Jl, 9m = 2 W • r • !; ' 2 

are not in involution with respect to the canonical Poisson bracket defined by (2.1) [11]. 
2.5 Rational Ruijsenaars— Schneider model 

Let us consider tensor field A, which is proportional to canonical bivector P 

A = (qx H h g„)P, p = qi H h g n . 

In this case equations (2.7) have the following partial solutions 

J k = ytxL\ fc = ±l,±2,...,±n, (2.16) 
where L is the Lax matrix of the Ruijsenaars-Schneider model 

n / 2 \ i/ 2 

As above recursion operator produces only the Hamilton function. It is easy to prove that traces 
of powers of the Lax matrix L (2.16) satisfy to the following relations 

PdJ± 2 = - r P'd In J k , k = ±l,...,±n, (2.17) 
k 

instead of the standard Lenard-Magri relations [25, 30]. Moreover, similar to the Calogero- 
Moser system, there are other solutions K m of these equations (2.17), which are described in [4]. 
Remind, that the so-called principal Ruijsenaars-Schneider Hamiltonian has the form 

.,2 \ 1/2 



1 / 2 \ 

jr »-i ( * +J -'-B™ ,, ** ) S( 1 "te^?) 



k=l j+k 

and that the rational Ruijsenaars-Schneider system is in duality with the corresponding variant 
of the trigonometric Sutherland system, see [4] and references therein. 

We want to highlight that for all integrable systems listed in [24, 30, 31, 32] the second 
Poisson bivector P' (1.1) is a Lie derivative of the canonical Poisson bivector P along the vector 
field Y = AdH (1.2), where tensor field A usually has a very simple form. 

In the next section we show that such simple tensor fields A may be useful to search for new 
integrable systems. 



8 



A.V. Tsiganov 



3 Noncommutative integrable systems 

The extreme rarity of integrable dynamical systems makes the quest for them all the more 
exciting. We want to apply tensor fields A to partial solution of this problem. Below we present 
a method to construct a new family of three dimensional noncommutative integrable systems. 
Let us consider natural Hamilton function on M = M 2n 

n 
i=l 

and bivector A associated with the rational Calogero-Moser system (2.12) 

A = (piqi H \-p n q n )P, 

where P is canonical Poisson bivector (2.1), (2.5). 

In previous section we have discussed partial solutions of the equations (2.7), here we want 
to discuss their complete solution. 

Proposition 1. The Lie derivative of P (2.1) along the vector field Y 

P' = C Y P, Y = (p igi + ■■■+ Pn q n )PdH (3.1) 
is a Poisson bivector compatible with P if and only if 

H = £,$ + \f(»,K..M. (3-2) 

~ q{ \ii qi qij 

Here F is an arbitrary homogeneous function of zero degree function depending on the homoge- 
neous coordinates 

_<12 _ <fy _ qn 

X\ — , X2 — , • ■ ■ , X n — 1 — 

qi qi qi 

The definition of the homogeneous coordinates may be found in [13]. Proof is a straightfor- 
ward calculation of the Schouten bracket (2.7). 

It is easy to see that some Hamilton functions separable in spherical coordinates and Hamilton 
functions for the rational Calogero-Moser systems associated with the A n , B n , C n and D n root 
systems have the form (3.2). 

We got accustomed to believing that the notion of two compatible Poisson structures P 
and P' allows us to get the appropriate integrable systems [15, 25, 30, 31, 32]. In our case 
recursion operator N = P'P^ 1 reproduces only the Hamilton function 

triV fc = 2{2H) k . 

It allows us to identify our phase space M = R 2n with the irregular bi-Hamiltonian mani- 
fold [25, 31], but simultaneously it makes the use of standard constructions of the integrals of 
motion impossible. 

We do not claim that all the Hamilton functions (3.2) are integrable because we do not have 
an explicit construction of the necessary number of integrals of motion. Nevertheless, even in 
generic case there is one additional integral of motion. 

Proposition 2. The following second order polynomial in momenta 

C = {piqi + ■ ■ ■+p n qn) 2 -(«? + ••• + ql)H 

is a Casimir function of P' , i.e. P'dC = 0. 
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Consequently we have 
{H,C} = 0. 

It is enough for integrability at n = 2 when we get Hamilton functions 

separable in polar coordinates on the plane. 

At n > 3 we can make some assumptions on the form of the additional integrals of motion. 
For instance, let us postulate that our dynamical system is invariant with respect to translations, 
i.e. that there is a linear in momenta integral of motion 

flpost = Pi H \-Pn, {H, impost} = 0. 

It leads to the additional restriction on the form of the proper Hamilton functions (3.2) 

1 ( 93 - 92 94 - 93 9n - 9n-l 



_ i yl2 ~ 9iJ V92 - 91 92 - 91 91 - 92 

which generate bi-Hamiltonian vector fields 

X = PdH = P'd In iJp^ t (3.3) 

equipped with the four integrals of motion 

H\ = impost; H 2 = H , H3 = C, i?4 = {Hi, C} (3.4) 

with the linearly independent differentials c?H. According to the Euler-Jacobi theorem [19] it 
is enough for integrability by quadratures at n = 3. 

Remind that the Euler-Jacobi theorem [19] states that a system of iV differential equations 

±i = Xi(xi, . . . ,x N ), i = l,...,N, (3.5) 

possessing the last Jacobi multiplier [i (invariant measure) and iV— 2 independent first integrals is 
integrable by quadratures. In our case N = 6, we have four independent integrals of motion (3.4) 
and fj, = 1. 

So, at n = 3 the following Hamilton functions 

H 2 = pf + pi + p| + 7 — r^-G f *^* > ) (3.6) 



(92 - 9i) 2 V92 - 91 

labelled by functions G generate integrable by quadratures Hamiltonian equations of mo- 
tion (3.3)-(3.5). Because 

{Hi, H 2 } = 0, {H 1} H 3 } = H 4 , {Hi, H 4 } = 2Hf - 6H 2 , 

{H 2 ,H 3 } = 0, {H 2 ,H 4 } = 0, {H A ,H 3 } = 4H 1 H 3 (3.7) 

we have noncommutative integrable systems with respect to the canonical Poisson bracket, see, 
for instance, [21] and references therein. 

Of course, in the center of momentum frame, the total linear momentum of the system is zero 
Hi = and we have three integrals of motion H 2 , H3 and H4 in the involution that is enough 
for integrability at n = 3 and n = 4. 



10 



A.V. Tsiganov 



On the other hand, Hamilton functions (3.6) define superintegrable systems in the Liouville 
sense 

= 0, i,j = 1, ... ,4, 

with respect to the second Poisson bracket {•, •}' associated with the Poisson tensor P' (3.1). If 
we put 

we can obtain a well-known Hamiltonian for the rational Calogero-Moser system (2.13) 

rr 2 , a 2 . g 2 g 2 

-«CM = > Pi + 7 To + 7 + 7 ^2" 

In this case there are other polynomial integrals of motion (2.14) and other Casimir functions 
of P' (2.15). This system was separated by Calogero [7] in cylindrical coordinates in M 3 . Being 
a superintegrable system it is actually separable in four other types of coordinate systems [5]. 
These variables of separation may be easily found using either the generalised Bertrand-Darboux 
theorem [35] or methods of the bi-Hamiltonian geometry [16]. Remind, that variables of sepa- 
ration are eigenvalues of the Killing tensor K satisfying equation 

KdV = 0, (3.8) 

where V is potential part of the Hamiltonian H. 

Any additive deformation of this function G(x) leads to the integrable additive deformation 
of the rational Calogero-Moser system, for instance, if 

G(x) = G{x)+ a , 
x 

then one gets an integrable system with the three-particle interaction 

~ a 
Hcm = Hcm + 



(<?1 - Q2) (<72 - 93)' 

For this Hamilton function we couldn't find any polynomial in momenta integrals of motion ex- 
cept Hi, H3 and H4 (3.4)). Moreover, we couldn't get variables of separation using the standard 
(regular) methods such as generalised Bertrand-Darboux theorem [35] and bi-Hamiltonian al- 
gorithm discussed in [16]. Namely, in contrast with the case a = at a 7^ the Killing tensor K 
satisfying (3.8) has only functionally dependent eigenvalues. 

At n = 3 in order to get rational Calogero-Moser systems associated with other classical 
root systems and their deformations we can postulate an existence of the fourth order integral 
of motion 

n 

# P ost = J2p 2 p 2 + Mq)pt + g(q), 

ij^j k 

with some unknown functions fk{q) and g(q)- However we do not have an exhaustive clas- 
sification as of yet. 

In generic case at n > 3 we can use other hypotheses about additional integrals of motion 
commuting with H (3.2). 
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Of course, construction of such integrable systems is trivial and closely related with con- 
struction of the group invariant solutions of partial differential equations through imposing side 
conditions [27, 28]. Remind, we can look for solution W{q) of the Hamilton-Jacobi equation 

2 



H (P^)=J2[^q i ) +W( qi ,...,q n )=S, Pj 



dW 



dqj 

up to the side condition 
S(q,p)=0. 

If {H, 5} = f(q,p)S then this side condition is consistent with H and the corresponding integrals 
of motion are defined modulo 5 = 0, i.e. 

{H,H k } = g k (q,p)S. 

Here / and g k are some functions on phase space and W{q) is the so-called characteristic 
Hamilton function. 

In our case the side condition is related with the transition to the center of momentum frame 
5 = i?p OS t = pi + P2 + P3 = 0, 

which is always consistent with the Hamilton function (3.6) and we have three integrals of 
motion H2, H3 and H4 (3.4) in involution by modulo 5 = (3.7). Construction of the variables 
of separation for the Hamilton-Jacobi equation with a side conditions is discussed in [6]. 
In quantum case we can consider the Schrodinger equation 

= H = A + V( qi ,...,q n ), 

where A is the Laplace-Beltrami operator on M = M 2n , and study solution of this equation that 
also satisfies a side condition 

5f = 0. 

The consistency condition for the existence of nontrivial solutions ^ is a standard 
[H,S] = fS. 

In this case linear differential operator K will be a symmetry operator for H modulo S^> = if 
[H,K]=gS. 

Here / and g are some linear partial differential operators, see examples and discussion in [6]. 

We assume that the quantum counterpart of H (3.6) could be embedded in this generic 
scheme. 



4 Conclusion 

We have demonstrated that the trivial deformations of the canonical Poisson bracket associated 
with the well-known integrable systems have a very simple form defined by some 2-tensor field 
A acting on the differential of the Hamilton function. We have shown a collection of examples 
and also proven that such tensor fields may be useful for searching new integrable by quadratures 
dynamical systems. For example, we have proven the noncommutative integrability of a new 
generalisation of the rational Calogero-Moser system with three particle interaction. 

In fact, we propose a new form for the old content and believe that this unification is a next 
step in creating the invariant and rigorous geometric theory of integrable systems on regular 
and irregular bi-Hamiltonian manifolds. 
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